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Abstract 

We describe heavy baryons as bound states of a quark and a diquark. For this purpose we 
derive the Faddeev equation for baryons containing a single heavy quark from a Nambu- 
Jona-Lasinio type of model which is appropriately extended to include also heavy quarks. 
The latter are treated in the heavy mass limit. The heavy baryon Faddeev equation is then 
solved using a static approximation for the exchanged quark. 
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1 Introduction 



The physics of hadrons with one heavy quark coupled to some hght degrees of freedom becomes 
much simpler in the heavy quark limit due to the arising heavy flavor and spin symmetries. 
The basic concepts have been worked out within the past few years 0] (see also p| for recent 
reviews). The heavy quark symmetries are realized approximately as long as the heavy quark 
mass niQ is much larger than a typical QCD scale, uiq 3> Aqcd, which is the case for Q = b,c. 
Phenomenologically, they result in an approximate degeneracy of heavy hadrons differing only 
in the heavy quark spin, and in a definite scaling behavior of observables with the heavy quark 
masses. 

Since the heavy quark inside the hadron has to be almost on-shell, it is useful to reformulate the 
heavy quark dynamics in terms of the off-shell momentum k = p — rnqv, where p is the total 
momentum and v the four- velocity (y"^ = 1) of the heavy quark, with k'^ -C rriQ. Furthermore the 
heavy quark spinor is split into large and small components by means of the projection operators 
and In the heavy mass limit {niQ oo) only the large component 

Qv{x) = l±ie*-'5--Q(x) , {^Q, = Q„) (1) 
survives, and the heavy quark part of the QCD Lagrangian reads 

CuQET = QSv-V)Q, + Oil/ mo) , (2) 

where is the covariant derivative of the strong interaction. Note that the heavy quark's velocity 
is a conserved quantity in the infinite mass limit. The operator [iv-T>) is obviously insensitive 
to spin and flavor quantum numbers of the heavy quark, therefore the corresponding hadron 
properties are determined by the light quarks only. The systematic expansion of QCD in powers 
of l/riiQ results in the so-called Heavy Quark Effective Theory (HQET) 0, ^. Recently, the ideas 
of HQET have been successfully extended to the sector with two heavy quarks (Heavy Quarkonium 
Effective Theory ) . Heavy fermion effective field theory techniques are also employed to derive a 
consistent expansion scheme for baryon chiral perturbation theory with SU{3)f flavor symmetry 
i- 

Although HQET provides the tool for calculating the perturbative QCD corrections to short- 
distance processes, one is left with the non-perturbative features of QCD that lead to the formation 
of hadrons (containing also light quarks) at low energies. As is well known, the low energy 
properties of light quark systems are dictated by the approximate chiral symmetry of QCD and 
its spontaneous breaking. The interplay of heavy quark symmetries and chiral symmetry has 
been formulated in terms of effective lagrangians for heavy meson fields 0, ^. An effective 
theory describing the low-energy interactions of heavy mesons and heavy baryons with Goldstone 
bosons TTjK,'!] was discussed in ref. |Ty]. Their results were applied to strong and semileptonic 
decays using spin wave functions of a non-relativistic quark model. The construction of chiral 



perturbation theory for heavy hadrons is also discussed in ref. where the chiral logarithmic 
corrections to meson and baryon Isgur-Wise functions are given. 



In ref. |T2| (see also |13| , |TJ]) we have presented an extension of the NJL quark model, consistent 
with chiral symmetry for the light quarks and the heavy quark symmetries, from which we derived 
an effective heavy meson lagrangian and predicted its parameters. In the present paper we shall 
use this model for the description of heavy baryons. 

The formation of baryons as bound states of quarks can be studied in general in two comple- 
mentary pictures. One is based upon the limit of infinite number of colors [Nq) in QCD, where 
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baryons arise as solitons of the meson fields |15|. In fact, heavy baryons have been recently |]T6| 



described as a bound state of a chiral soliton and a heavy meson within the extended (bosonized) 
NJL model of [jl2|. This so-called bound state approach [T^ has also been applied to a descrip- 
tion of heavy baryons within the Skyrme model |T^. In ref. []19|, the large Nq limit has been 
analyzed and an induced algebra with respect to the spin-flavor symmetry was derived. This 
model-independent approach was illustrated in the non-relativistic quark model and the static 



SU{2) chiral soliton model. Furthermore, in ref. pO[ a new formalism for treating baryons in the 
1/Nc expansion, based on an analysis of quark-level diagrams, has been given. 

On the other hand, for finite number of colors {Nc = 3), two of the three quarks in a baryon 
are conveniently considered as diquarks. Then baryons can be described as bound states of quarks 
and diquarks, resulting in Faddeev type of equations |2^. Here diquarks play a similar 

role as the constituent quarks, i.e. they serve as effective (colored) degrees of freedom for the low 
energy dynamics inside the hadron. 

In ref. the phenomenologically successful Nambu-Jona-Lasinio model has been converted 
into an effective hadron theory, where meson and diquark fields are built from correlated qq and 
qq pairs, respectively, and baryon fields are constructed as bound quark-diquark states. This 
approach, which is fully Lorentz covariant, also allows to study the connection between the soliton 
and the quark-diquark picture of baryons as well as the synthesis of both pictures 



In this paper we will apply the above mentioned NJL model with heavy quarks to a quark- 
diquark description of heavy baryons in lowest order of the (l/mg) expansion. For this aim we 



will generalize the hadronization approach of ref. |^ to the lowest -lying baryon states containing 
a single (infinitely) heavy quark, given by the following fiavor multiplets p5| , p6| , p7| : 

i) Three spin-1/2 baryons where the light quark spins are coupled to zero and which transform as 
3f with respect to the light fiavor group SU{3)f, represented by the antisymmetric matrix: 






>,/3=l/2 



Aq 




-1/2 




(3) 



Here is the (conserved) velocity of the heavy baryon, and (1 + '^)/2 is a projector on the 
large components of the baryon spinor. The two spin orientations of the spinors form the 
(in this case trivial) spin-symmetry partners. 

ii) Six spin multiplets containing a spin-1/2 baryon By and a spin-3/2 Rarita-Schwinger field 
B*^ transforming as 6^ under SU{3)f: 



v^B^ + B*/ , 



1 + 







1 

V2 



/3 = -l 



,^3=0 

, ^'(*) ^'(*) 

V "Q,/3=l/2 "Q,/3=-l/2 



"Q,/3=l/2 ' 
"Q,/3=-l/2 



(4) 
(5) 



Heavy fiavor symmetry leads then to a degeneracy of the residual masses AM = M — mq for 

different heavy flavors. 

The paper is organized as follows: 

In section Q we deflne our model which follows from the extended Nambu-Jona-Lasinio model of 
ref. by keeping those parts relevant for the formation of heavy baryons. Furthermore, using 
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the above quoted hadronization approach our model is transformed into an effective theory for 
heavy baryons. The resulting Faddeev equation is analyzed in section where we also introduce 
some approximations which facilitates its numerical solution. Numerical results for the heavy 
baryon masses are presented in section |[ Some technical details are relegated to the appendices. 



2 NJL model with heavy quarks 

2.1 Quark lagrangian and generating functional 

Let us consider a QCD-motivated NJL model with color-octet current-current interaction: 



^NJL = -3 G , (6) 



where = ^'7^^\E' , and = {q, Q) contains both the light q = u,d,s and the heavy quarks 
Q = c, b, and G is an effective coupling constant of dimension mass~^. This interaction can be 
Fierz-rearranged so that it acts only in the attractive color singlet quark-antiquark and color 
anti-triplet quark-quark channels [^. These channels can be considered as the physical ones 



in the following sense: Mesons are built up from color singlet quark-antiquark pairs, while in a 
baryon two of the three quarks form a color anti-triplet state. 

The interaction is chirally invariant. For the light quark flavors, the quark-antiquark interac- 
tion leads to a spontaneous breaking of chiral symmetry, which is accompanied by a dynamical 
generation of a constituent quark mass |2^. Further contributions from the mesonic sector will 
not be taken into account for the calculation of heavy baryons. Our effective (constituent) quark 
lagrangian is then defined by: 

£ = q{t^-m)q + Q,{tvd)Q, + , (7) 

where represents the diquark correlations arising from the Fierz transformation of -CJvjl ®- 
Since we are interested in baryons with a single heavy quark we will keep only the light-light and 
light-heavy diquark correlations, i.e. we ignore diquark correlations with two heavy quarks. Then 
the diquark interaction^ reads: 

C^' = 2G'i(rr°g)(gr>^) + (^ rr Q.) (0. g^) . 

Here, q'^ and q'^ are the charge conjugated spinors, given in the Dirac representation by 

q' = Gf, r = q'^G, C = 27270 = -C^ 

Furthermore, the vertices F of the diquark currents are defined by: 

(light diquarks), 
(heavy diquarks). 




2 



01 := {1, «75, ^{lt^-l^Vf,), - 'j^v^)-f5} , (8) 



^ Heavy quark symmetry and chiral symmetry allow for different coupling constants Gi, G3 for the light and 
heavy diquark sector. A factor 2 is absorbed into the vertices F^. 
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where ie^ denotes the color Clebsch-Gordan coefficients (given by the Levi-Civita tensor {€^)ij = 
^Aij) couphng the product representation 3 x 3 to 3*. Moreover, = a = 1, ... ,8 are the 

generators of the hght flavors SU{3)f group, and = | \/f The symbol denotes the Dirac 
matrices, where in the heavy sector the diquark degrees of freedom have been rearranged to make 
vector and axial vector diquarks transversal with respect to the heavy quark velocity v'^. 
Due to the Pauli principle , the light-light diquarks in the (symmetric) Gp flavor representation 
exist only as Lorentz axial-vectors, i.e. with = ij^/ ^/2, whereas the diquarks in the (antisym- 
metric) 3f representation occur with = {1, i'j^, ^7/^75}. For the vertices of the heavy-light 
diquarks^ there is no restriction since heavy and hght quarks have to be treated as different 
particles. 



2.2 Hadronization 

Our aim is to convert the quark theory defined by (0) into an effective theory of heavy baryons. 



For this purpose we follow the hadronization approach of ||22[ , where baryon fields were built from 
diquark and quark fields. 

Baryons with a heavy quark can be constructed from either a light diquark and a heavy quark 
or a heavy diquark and a light quark. Therefore we introduce both types of diquark fields into 
the generating functional 



Z = J Dq Dq DQy DQ^ exp 'i' J ^ 



with the help of the following identities]^: 



1 = J 1[Dk''Dk^''6{k^" -qr"q^)6{K" -q^r^q) 



g-r-g^) 



h.c. 



n Dk'' Dk^" DA" DA^" exp ; 

a 

I n Dk'^" Dnr" Sinr - Q, rr q') 5(«:."'' - ^ Qv) 



n ^4"" ^^r" exp - / [Ar {kC - q') + h.c. 



(9) 



(10) 



Note that the heavy diquark field A^, = (A^)^. (and (kJ,)^) with j = u,d,s transforms as a triplet 
under SU{3)f- Analogously, we also introduce two different baryon amplitudes for the heavy 
baryons: A baryon field amphtude built up from a light diquark and a heavy quark (A'„ ~ 
(g'^g) Qy) and an amplitude constructed from a heavy diquark and a light quark (y^ ~ {.<fQv) l)- 
These fields are introduced again by inserting the following identity into the generating functional: 



n DX'^ DX^ DX^ DX^ 



a 



^In the following we will refer to the diquark fields simply as light and heavy diquarks, respectively. 

•^We have defined the conjugated fields in such a way that (AT 70)^ := A'''r7o. Since some of the are 
anti-hermitian, the respective diquark fields behave likewise, leading to an unusual phase convention. As usual 
irrelevant normalization factors in front of functional integrals are omitted. The numerical factors 1/2, 46*1, 2G3 
in eqs. (|9|)-([T^) have been introduced for later notational convenience. In order to keep notations transparent, we 
sometimes write / instead of / d^x etc. 
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exp I 



[X:{x,y)-AG^Q,{x)[q{y)T^q\y)]] X:{x,y) + h.c] , 



/ n DY,^ 

■ exp i 



"° DY,-^ ny^t'' Dy,"'^ 



y,t^ix,y) - 2Gsq.{x)[QM^7 Qjiy)]] Y,t^{x,y) + h.c.} . (12) 



Note also that -^X^ = and f^3^^ = 0. For later convenience, we have introduced the baryon fields 
in two different flavor bases: The flavor index a of (contained in a = {a, A, a}), a = 0, . . . , 8, 
refers to the generators of the S'f/(3)i? group and the unit matrix, whereas the flavor indices i,j 



of X 



(ao 



{y4,a}), z,j = 1..3, denote the light flavor of the quark and the heavy diquark. 



respectively. 

Owing to the 5-functionals in (pD-(p!2D, we can replace the quark bilinears in the interaction 
d^) by the flelds k, k^, and n\ k\j, resulting in a lagrangian which is bilinear in the quark flelds. 
After integrating out the quark fields as well as the auxiliary fields k, we obtain an effective theory 
in terms of the diquark fields A, A„ and the baryon fields X^, X^.Y^.y^ (see ( |A.1| )). This effective 
baryon theory is highly non-local and we will work it out in the low-energy regime. Since we are 
interested in the description of individual baryons we expand the corresponding effective action 
up to second order in the baryon source fields^ X'^^Y"^. This yields (see appendix A for more 
details): 



Z 



J DA DAt DA„ DAl DX^ DX^ ■■■ exp {iS,s) , 



Tr log g;^ - i Tr log g-' - - Tr log(l - g A~g A^^) 



{G 
{G 



B XX 



B YX 



['^B XY) i'f 

fry \l3ol3o' 
[^B YYjm'-ii 



V 



For this purpose we have defined the inverse of the light and heavy quark propagators as: 



9 ^ix,y) 



{%^, - m) b{x -y)^ \ \C Al{y) g,{y, x) A,{x) C] 
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C^g^C 



1 + 



'w-d^)6{x -y) 



9v 



C'9lC 



(13) 

(14) 
(15) 



with the abbreviations A := A" F", At := At"F", A^ := A^"» F^o, Aj := Aj"" F"o. The symbol 
Tr denotes the trace over coordinate space and fiavor, color and Dirac indices. The entries for 
the matrix Gb are of 0*'^ order in the baryonic fields Xy and Y^, but still contain all orders of 
quark-diquark contributions, see ([A.7| ), ( |A.8| ). After the Gaussian integration over the auxiliary 
baryonic fields X^ and Yy, the matrix Gb can be identified as the propagator for the baryon fields 
X^ and y^. 



2.3 Derivation of the baryon propagator 

To obtain the effective theory of the heavy baryon fields, we have to integrate out the diquark fields 
A. This can only be done in an approximate fashion due to the presence of the quark loops coupling 

^ Higher powers of baryon fields in (|l^) would lead to the generation of baryon interactions in the effective 
lagrangian which are outside the scope of this paper. 
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to the diquark fields in the effective action (plsl) . First we derive the free diquark propagators 

in leading order of the loop expansion, i.e. we expand the term Tr log [1 — g ^ g A"!") in (|1^) up 

to second order in the diquark fields. This yields: 

D^^'-\x,y) = --Ls^^' 6{x-y) + '-tT\T^ go{x,y)T''' Uy,x)] (16) 

{D,)ff^\x,y) = -±-6^-^-' 6,y6{x-y) + '-tT\r^,'^ g,{x,y)r^/ {~go),Ay^^)] (17) 



for light and heavy diquarks, respectively, with {go)ii'{x, y) := {i(^x — i^iY^ S{x—y) Suf denoting the 
free light quark propagator. The diquark polarization tensors (the terms in the square brackets) 
are calculated in appendix P and graphically depicted in Fig. [l|. Note that the heavy diquark 
propagator is a diagonal flavor matrix, {Dy}^^, = diag(D", -D^, -D^). 





Figure 1: Diquark polarization tensors for light and heavy diquarks, generated by the loop expan- 
sion of the quark determinant. 

Now the functional integration over A and A^ will be performed in leading order of the cluster 
expansion 
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(CO 1 \ 
E ((n)j = exp {(F)a + 0{{F')a - {F)l)} , 

where the functional average over the diquark fields is defined by: 

/ DAW AD AID A, [. . .] exp i J (At"(D-i)"'3A/3 + Ar°(/^,"')°°*A(^o 
^^■■■^^"^ ~ J DAWADAIdA, exp i/(At"(Z}-i)"/3A/3 + Ar"(i^-i)°o^oA(^°) 

Obviously, the diquark propagators ([T6|) and (|1^) are then given by: 

zD^^\x,y) = {A^{x)A^^\y)) , (19) 
z{D^)ff{x,y) = (A4«(x)At*'(y)) . (20) 

Now we are left with the task to evaluate the functional average (Gb) over the diquark fields. 
With the help of Wick's theorem this functional average can be reformulated into the sum over 
all possible contractions of the diquark fields in into free diquark propagators (AA'''). Hereby, 
diagrams with an arbitrary number of intermediate diquark propagators arise. This refiects the 
fact that even in the constiuent quark model the baryon is in principle a many-particle system, as 
soon as the Dirac sea is included. In the following, we will only keep the simply nested exchange 
diagrams shown in Fig. 13. Here the first diagram describes the independent propagation of a quark 
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(Gb) 



+ 





Figure 2: Diagrammatic expansion of the baryon propagator. 



a 1 [ 



b) n 





d) 




Figure 3: a) The unperturbed baryon propagator describing the independent propagation of a 
quark and a diquark. By construction of the diquark propagator, this diagram includes already 
the diagram shown in b). c) Exchange diagram to b). d) The exchange diagram c) redrawn to 
show the equivalence to Fig. ^. 

and a diquark. A consistent treatment requires to include also the corresponding exchange diagram 
shown in Fig. Notice that there arise also self-energy corrections for the quark propagator due to 
inserted diquark propagators which can be absorbed into the constituent quark mass. Moreover, 
there emerge higher order contributions to the exchange diagram in Fig. |^c) due to corrections to 
the quark-diquark vertex (see Fig. Ha). A consistent treatment of the baryon would then require 



a) 




b) 




Figure 4: a) A typical higher-order correction to the quark-diquark vertex, b) Corresponding 
higher-order correction to the diquark propagator, which is discarded in this paper. 

to include the same correction into the diquark propagator. This complicated procedure which 
goes beyond the one-loop expression for the diquark propagator is clearly outside the scope of this 
paper. 

Let us remark that in the case of light quarks the class of diagrams shown in Fig. □ just 
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represents the minimal subset of diagrams which is required to fulfill the Pauli-principle in the 
quark-diquark picture of baryons p2| , [2^ . The above series of diagrams defines a kind of ladder 
approximation which in the next step will be summed up to yield a Faddeev type of 3-body 
amplitude. Hereby, the presence of both light and heavy diquarks induces a 2x2 matrix structure. 
Therefore, we introduce the matrix of the free quark-diquark propagation G = diag {Gxx, Gyy), 

{Gyy)^]%,{x,x';y,y') = imf/^^x' ,y'){g,),e{x,y) . (21) 

Here, and gQ denote the free propagators of heavy and light constituent quarks, respectively. 
Now we define the interaction matrix if, whose elements mediate the transitions ^ and 
yv ^ yv As can be read off from Fig. |^ by straightening the diquark lines, these transitions are 
mediated by quark-exchange. In Fig. ^ the elements of H are depicted graphically. We have: 



H 



{Hxy)%{x, x'- y, y') = - T^/ (^o(x, y)r^),v 5{x' - y) 5{y' - x) 



HxY 
Hyx Hyy 

1 
2 

{HYx)^f\x,x';y,y') = 1 (r^'^o(x, y))., - 5(a; - y') , 

{HYY)^f/Ax, x'- y, y') = ^ t^/ y) 5{x' - y) 5{x - y') 5,, • (22) 




A. A 5 A„ ? 

HxY Hyx Hyy 

Figure 5: Elements of the interaction matrix H . 

Here, F^, = F^ C . Note that Hyy formally mediates the exchange of the heavy quark. 
With the above definitions of G and if, the baryon propagator defined by the series of diagrams 
shown in Fig. ^ becomes: 

{Gb) = G Y.{H Gf = G {I- H GY^ . (23) 

k 

It has the standard form of a Faddeev type of 3-particle propagator. 
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3 The Faddeev equation 



Baryon masses are given by the poles of the baryon propagator (PBD, which results in a Faddeev 
type of equation for the baryon wave functions: 



1-GH) =0. (24) 




p/2-q 




p/2 + q 



p/2 - q 



y "0 




p/2 + q 



p/2 + k 




p/2 + k 



p/2-k 



n-{k + q) + 



p/2 + q 




p/2-q 



ii-(A; + g) 



p/2 - k 



p/2 + q 



Figure 6: Graphical representation of the Faddeev equation for heavy baryons. 

For the solution of this equation it is convenient to switch to momentum space. Denoting the 
total momentum of the baryon (with rnqv subtracted according to the phase factor in (|l|)) by p 
and the relative momentum by 2q, 2k, respectively, we obtain (see Fig. 13): 



X:{p,q) 



2 Jk 

- I 

2 Jk 



+ - 



D^'^'ip.) g.ip+) + q) T"'),. 3^4°(p, k) , (25) 

(D,)^,"°'(p_) ((7o)n'(p+) (r^ Uk + q))vr rr' K^iv.k) 

(DXr'iP-) i9oh'{p+) f ^° ~9v{k + q) f 3^4v(P> k) . (26) 



Here p+ and p_ are abbreviations for (p/2 + q) and {p/2 — g), respectively, and stands for 
/ (i^A;/(27r)^. Note that in momentum space one has g^i^k + q) = C^gQ{ — {k + q))C (here, the 
transpose refers to Dirac indices only). In the following, we will restrict ourselves to the lowest 
lying diquark states, namely the scalar and axial-vector diquark channel. 

Furthermore, we will not perform a full dynamical calculation of the Faddeev amplitudes, which 
would require a tremendous amount of numerical work. Instead, we will consider the diquarks 



as point-like bosons. This reduces the Faddeev equation (|25|) , (|26D to a Bethe-Salpeter type of 
equation. For the light diquarks, we use Klein-Gordon propagators D{k'^) = Z[k^ — M^]~^, with 
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masses M = M^, and wave function renormalization factors Z = Z^^ Za-, where the indices refer 
to scalar and axial-vector diquarks, respectively. 

Accordingly, for the heavy diquarks we use the propagator of a free bosonic field which has 
been subject to the heavy quark transformation i^: 

[D^^'tUv k) - J , X ( ^ for scalar diquarks 

vk-AMj+ie " Sf/i^j' - f V for axial-vector diquarks 

Here, AM = M—uiq is the residual diquark mass and Z denotes the wave function renormalization 
factor. The heavy scalar and axial vector diquarks are degenerate due to the heavy quark spin 
symmetry. The diquark masses and Z-factors appearing in the above expressions are calculated 
using eqs. (p!6| ) and (|1^) (see appendix ^ for details). Numerical values are given in section ^. 

The color structure of the interaction matrices H can conveniently be rewritten in terms of 
the color projection operators onto singlet and octet states: 

ttAB _ '^^CD ^^DE _ _{l)pAB , \{fi)pAB W pAB _ 

a/2 ~ 2 ' ~ 3 ' 



Then eqs. (|25|) , ( ^61) can be shown to decouple in color space and the parts describing the relevant 
color singlet amplitudes can be identified. 



In the case of exact flavor symmetry, the flavor structure of eqs. (^51) , (^) is trivial, since all prop- 
agators are proportional to the unit matrix in the (light) flavor space. Then the flavor matrices 
connected with the fields X°- project out either the symmetric {Qp) or the anti-symmetric (3^?) 
flavor part of 3^jj, respectively, which leads to a total decoupling oi Qp and 3i? amplitudes. 
The explicit breaking of the SU (3) p flavor symmetry induces a mixing of these amplitudes. As a 
first approximation we will keep the flavor structure in the kernel of eqs. (p5l), (^Bf ), but neglect 
the mixing effects for the amplitudes. The error involved is of first order in (m^ — m^) for the 
baryon wave functions, i.e. of second order for the mass spectrum. 

Furthermore, we are interested only in ground state baryons, which are given by s-wave (quark- 
diquark) amplitudes as discussed in For s-wave baryons, the transversal momentum of the 



exchanged quark should be of minor importance and will be neglected. This approximation implies 
the replacement: 

}^^{vk)ij Uk) ^ T^— -' (27) 

—v-k — m + te 

where the relations = and = ^ Ty were used. Now we insert ( ^7]) into eqs. (|25|) , 

([26| ) and perform the integration over the relative momentum q, defining X{p) = JqXy{p,q) and 

yvip) = Iq 3^1) (p, q)- This leads to a decoupling of the terms ^ and the longitudinal parts v^X^ 
(where X'^ denotes an amplitude containing an axial-vector diquark). Therefore, the equations 
now enforce the Bargmann-Wigner condition, i.e. to leading order in l/mq the baryon amplitudes 
are eigenstates of the projector and the transversality condition (^ is fulfilled (see for 
more details about these conditions). 

Finally, linear combinations of the baryon amplitudes have to be found which are eigenstates 
of projection operators onto spin 1/2 and spin 3/2. For this purpose, we introduce the following 
set of spin projection operators: 
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Here, ^^^'^^P^y projects onto spin 3/2 while t^t^ projects onto the part of spin 1/2 which is transver- 
sal with respect to v^. The Dirac structure in (p5D, (pGf ) stemming from the interaction vertices 
can now be rearranged as0 



The system of integral equations (pSj), (|26| ) can then be diagonalized with respect to the spin 
structure by introducing the following linear combinations: 



V 1 



[3/2] := ^ [=^/']p,^A'; . (30) 
v2 

and 

W2]y*f. := [3/2]pA^3;- . (31) 

Note that in the case of '^/2'3^t, and ^^^'^'^y'y both the heavy scalar and heavy axial vector diquarks 
contribute. 

Now eqs. (p5|) , (|26|) can be decomposed into two separate sets of coupled equations for the baryons 
with the light flavors in the 3^ representation (Ag, Eq) and the Qp representation (Sg'', 
fiq^), respectively. The equations for the fields ^^^'^^X*"', ^^^^^y^^j indeed turn out to be identical 
to those for the fields ^^/"^^X^", ^^^"^^y^^j which reflects the de generacy due to the heavy quark spin 
symmetry. 

The fact that the system is now diagonal with respect to the flavor and the spin content shows 
that the Pauli principle (which was inherent in the amplitude Xy containing the light diquark) is 
established dynamically in the amplitude y^. As shown above, the physical baryon which will be 
a linear combination of Xy and y^, always has a definite relation between flavor and spin quantum 
numbers. 



4 Numerical results 

Even with the simplifications introduced in the previous section, a complete solution of the re- 
sulting coupled integral equations for the baryon amplitudes Xy{p), y^ip) would be a rather com- 
plicated numerical task. 

To get a first rough estimate of the baryon spectrum, we apply an even stronger final approxima- 
tion which has proven to work quite well in the light baryon sector (compare 0, ^ 0)- This 
approximation implies the neglect of the total momentum dependence of the exchanged quark in 

^ Note that since all amplitudes are transversal now, the vertices can be brought into a similar form as 0^. 
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if, i.e. go{k) — 1/m with m being the hght quark mass, and is referred to as a static approxi- 
mation. Here, we will use the same approximation also for the propagator of the exchanged heavy 
quark, since, as can be seen from eq. (^), it differs from the light quark propagator only by a 
constant shift in the relative momentum of the amount mv'^. This shift can be absorbed into the 
definition of the wave functions entering eqs. (|25|) , (|26|) . 

With these approximations for the exchanged quarks, the integral equations (p5|) , (p6|) for the 
baryon amplitudes become a set of completely decoupled algebraic equations for the masses of spin 
1/2 baryons in the Sp representation and for those in the Gp representation, the latter including 
the spin 3/2 baryons. The final equations are of the form 



1 . 



(32) 



where Fij{vp) and = 1/2, 3/2, defined in appendix are integrals over the relative mo- 

mentum of the quark and diquark propagators. These integrals are UV-divergent. In a complete 
analysis of eqs. (|25|), p6|), the integration over the relative momentum q would be limited to the 
width of the wave-function appearing in the integral, which would render the integration UV-finite 
(compare e.g. [0). Here, we will employ an UV-regularization for the quark-diquark relative 
momentum integrals, where for simplicity the same method will be used as in the calculation of 
the diquark masses (see appendix ^). 

The usual NJL model parameters are fixed in the light pseudoscalar meson sector and the 
light baryon sector. For the latter we refer to the NJL model calculation in [^. In this work, 
a constituent mass of = nid = 450 MeV was used, from which the NJL cut-off and the 
strange quark constituent mass are calculated by a fit to the light meson spectrum, resulting in 
A = 630 MeV and rUg = 650 MeV. Moreover, there it was found that the coupling constant for the 
light axial vector diquark has to be increased relative to that of the scalar sector: Gia = 1.5 Gi. 
In the following, the value of the coupling constant in the heavy sector G3 is left equal to that of 
Gi- For convenience, the corresponding mesonic coupling constant Gmeson = 19.4 from the 
above fit is used ClS db normalization for Gi. 

First, some diquark masses for different values of Gi are presented in Table 0. The results for 
the masses of heavy diquarks are similar to those of heavy mesons |T^ . 
In Fig. ^ diquark masses and residual masses of the non-strange heavy baryons Ag and Sg 



G\/ Gijioson 












AM„ 


AM, 


1.1 


705 


895 


875 


1050 


1215 


360 


530 


1.2 


680 


870 


865 


1035 


1200 


345 


520 


1.3 


650 


845 


855 


1020 


1185 


330 


505 


1.4 


620 


820 


845 


1010 


1170 


320 


495 


1.5 


595 


795 


835 


1000 


1160 


305 


480 



Table 1: Light scalar and axial vector diquark masses and heavy diquark masses (all in MeV) for 
rriu = 450 MeV, m, = 650 MeV, A = 630 MeV, G3 = Gi, Gia = 1.5 Gi, G^^^on = 19.4 A-^. 

following from eq. (^) are shown as a function of the coupling constant Gi. The circles indicate the 
masses of the diquarks found in the fit to the light baryonsQ which determines the parameter 

^In this work a different regularization method was apphed, therefore the circles appear at shghtly distinct 
values of Gi. 
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Figure 7: Light diquark masses M^, M^, heavy diquark residual masses AM and heavy baryon 
residual masses AMaq, AMs^ as functions of the coupling strength. The mass splitting between 
Ac and Sc found in experiment is marked with an arrow. 

range to be used in the following as Gi/Gmeson = 1-2 — 1.3. Indeed, here we find a reasonable 
mass-splitting of 160-175 MeV for the heavy baryons Aq and Eg. The experimental value for the 
charmed baryons is M^^ — Ma^ = 170 MeV which is indicated with an arrow in Fig. |^. 



The following masses of charmed baryons have been determined so far |]3J, ^ 



Ac(2285) , Sc(2455) , 2^(2470) , ^]c(2705) , S*(2644). 
In the bottom sector only the data for A{,(5641) is reliable up to now.[| 

Our results in Table ^ were obtained without changing the parameter setQ used for the calculation 



of light meson and baryon spectra in |31]. For a value of G\ = 1.3 Gmeson, the mass splittings 



Ac and Sc — Ac obviously agree quite well with the experimental data. Note that the (pre- 



liminary) result Sfe — Afo = 173 MeV I^G] equals the corresponding mass splitting in the charmed 
sector. The Qc — Ac splitting comes out somewhat too large, whereas the splitting S* — Sc is 
somewhat too small. 

However, our calculation without l/mq corrections cannot account precisely for the mass split- 
tings within the multiplets. Using spin-weighted averages for these multiplets (with prelimi- 
nary results or estimates for the candidates not yet determined) would significantly improve 
our predictions for and (whereas the prediction for {S', S*} would probably be 

somewhat too small then). But a definite analysis requires more reliable experimental data, and 
perhaps even a re-analysis of the so far performed classification of heavy baryon states 



^There is first experimental evidence for Sb(5814) and E J (5870) 
^For comparison, also results for Gi = 1.4 Gmeson are shown. 
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^-^Xj meson 


So - An 


Sn — An 


VLn — An 




1.20 


162 


192 


495 


136 


1.30 


175 


198 


513 


146 


1.40 


189 


204 


531 


156 


Exp.: 


170 


185 


420 


174 



Table 2: Heavy baryon mass splittings in MeV. The NJL coupling constant Gi, which triggers 
the diquark properties, is given in units of the coupling constant of the (light) pseudoscalar meson 
sector. 

5 Summary and outlook 

In this work we have derived the Faddeev equations for heavy baryons in the quark-diquark 
picture, starting from a current-current interaction of the NJL-type model, where heavy quarks 
have been included in the heavy mass limit. 

After extensive use of functional integration techniques, the baryon propagator has been ob- 
tained as a series of quark exchange diagrams, which could be summed up to yield the bound 
state equations for the baryon wave functions. 

An essential feature of our approach is that two different types of quark-diquark configurations 
contribute to the baryon wave functions: One, where a heavy quark is coupled to a light diquark, 
and the other where a heavy diquark is coupled to a light quark. The properties of the heavy 
diquark are also fixed by the underlying NJL model and discussed in some detail, showing a similar 
mass pattern as in the case of heavy mesons. 

For a further analysis of the baryon bound state equations, we approximated the diquarks as 
point-like particles, which reduces the Faddeev equations to effective Bethe-Salpeter type of equa- 
tions. In an s-wave approximation (that means in our case to neglect the transversal momentum 
part of the interaction), the Bargmann-Wigner condition and the transversality condition for the 
heavy baryon wave functions are fulfilled exactly. Projecting onto color singlet states and defining 
the proper spin projections, we then observe a decoupling of the equations for suitably chosen 
linear combinations of wave functions. They describe the predicted multiplets of heavy quark spin 
symmetry: the lowest-lying spin 1/2 state {aSp multiplet in the light flavor's SU{3)f) and the 
spin 1/2 state (a 6i;' multiplet) which is degenerate with the spin 3/2 state. 

The Pauli principle for the two light quarks inherent in the light diquark system is established 
dynamically for the baryon wave functions: We have shown that in the SU{3)f limit baryonic 
states (with definite spin) where the light quarks are in the 3i?-representation do not mix with 
those states where the light quarks belong to the 6^ multiplet. This feature is preserved even if 
SU{3)f breaking is included in leading order. 

For a numerical estimate of the mass spectrum, a static approximation for the quark-diquark 
interaction was found, which treats both heavy and light quark exchange on an equal footing 
and leads to analytically solvable equations. By applying a parameter fit obtained from the light 
baryon spectrum, we then arrived at reasonable masses for the ground state heavy baryons. 



As we have shown in ref. |3^, the quark-diquark picture is also useful to estimate the couphng 
of heavy baryons to pions and weak heavy quark currents, defining the Isgur-Wise form factors. 
The results of the present work are to be extended in this direction in the future. 
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A Deriving the effective baryon lagrangian 



After inserting the functional constants (p|)-(p^) into the generating functional, the effective action 
reads: 



5, 



eff 



q [i (j) -m) q + Qv {iv ■ d) Qi 



(A.l) 



\j (g A" r" + At " r" g + g^F^ A^^g^ + (f F^AJ "°g, 

2 Jx,y ^ 



At"(x)A"(x) - / QMX"{y,x)A^"{x) 



+ 



8Gi Jy 

Jy' 

^ A7(a;)Aj,"«(x) - / q{y) f x)A|; 



+ 
+ 



x,y 



Ae°(x)F°(i/,x) r q{y) - 4G3 / q{y) r Y^{y,x)Y''{y',x) T q{y' 

Jy' 



xx + xx + Yy + yY. 



x,y 



The heavy quark fields are removed from the generating functional by a simple Gaussian integra- 
tion. For the terms containing light quark fields we employ the Nambu-Gorkov formula: 



/ Dq Dq exp i ]- j (g, g 

J Z, Jx,y 



A B 
C D 



(g' ) ^ llS'''' + 



exp 



iv, -V ) 



A B 



-1 



C D 



-T 
-T] 



(A.2) 



were in our case the components read: 



A{x,y) = ii^-m)S{x-y) - AGs Yaa,{x, z)Ya'a,{y, z) 



B{x,y) 
C{x,y) 
D{x,y) 

rj{x) 



-^[C Al{y)S,{y,x) A,{x)cf 

-A{x)5{x -y)C , 
-C A{x)6{x - y) , 
-A^{x,y) 



2e I A, 



X°'(x',2)A°'(z) S^{x',x) A^(x) C 



t Oto 
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2 / Ay(z)Y''(x,z) r 



C A,{x) S,{x,y) At °'(^)X"'(y,z) 



y,z 



(A.3) 



Here, Tr denotes the trace over all indices, including space-time indices. We define the inverse of 
the following matrix 



and the quantity 



A B 
C D 



1 + ^ 



gii 9i2 

921 922 



S;\x,y):=itvd)—^5{x-y) - 8G, X^{x,a) X"{y,a) . 

Z J a 

Note that the projection operator ^ is always present next to a or field. 



(A.4) 



(A.5) 



The effective action now reads: 



i S, 



cff 



8Gi 



I At "A" - ^ I Ar At -0 + Tr log S;' 

i 



+ -Tr lo 
2 



A B 
C D 



( 9ii 


9l2 \ 


( -V ) 


\ 921 


922 1 





+ 



X°(x,a)A"(a) S^{x,y) At °(a') X"(?/,a') 

/ xx + xx + Yy + yY. 

■J x.v 



x,y,a,a' 



(A.6) 



Expanding this expression up to second order in the baryon fields, we obtain the entries for the 
matrix Gb defined in ( pTSD as given by 

{Gb xxf^\x,x'; y,y') = A^x') g,{x,y) A^^\y') 

+^ J d^ad^z Af'ix') g,{x,z) A.^iz) C g22{z,a)kiC A^ia) g,{a,y) A^"' {y') , 



(Gb xYfv]'{x-,x'\ y,y') = ^ y A^{x) gvix,a) A^k{a) G g2i{a,y)kv Afjl' {y') 
{Gb Yxtf\x,x'- y,y') = \ f d'a A,f{x') g,2{x,a)uG Al{a) g.,{a,y) At^'(yO 



(Gb yy)5,e°;(x,a;'; y^y') = A^f{x') g,,{x,y)., Al^;\y') 



(A.7) 



where 



911 

922 
9l2 
921 



l-gA-gA^^r' 



9 , 



T 

-9ii 



{l-g Ag A^Y'g AgG^ , 
G\l-~gA^gA)-^~gA'^g. 



The expression g is defined as ^ = Ct g^ G, where {}^ means 'transposed in all indices' 
(including coordinate space). Note that the Dyson expansions for the terms gu yield quark- 
diquark contributions up to any order. 
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B Calculation of integrals in the proper— time scheme 

As far as the integrals (in momentum space) for the diquark propagators (|TB|) and (|T^ are con- 
cerned, the proper-time regularization procedure is apphed to the quark determinant appearing 
in (IT^) , such that Ward-Identities are respected. The essential techniques are outhned in [|^. For 
the integrals with one heavy and one light particle in the loop, the basic replacement reads (after 
Wick-rotation): 

ds e-^(« +™ ) — . / ds e-^(« +™ ) . (B.l) 



qj. _|_ 



We define 



(27r^) J [q^ ~ m? + ie) {v-q + v-p + ie) 
167r2 V ^ ^ 2' A2 ' 

+^ /'^rffl ^'-f-^'V (B.4) 

i r^s d'^q ( I q^ - 



(27r'*) J [q^ — m? + ie) {v-q + v-p + ie) \ 3 m 
TT {w? — iy-p)'^)^ ^ , 3 rn? — {v-p)' 



2 



327r2 m? '2' A2 ' 



Here V{a,x) = dt e denotes the incomplete gamma-function. 

B.l Diquark polarization tensors 



Starting from the expression (|TB|) for the inverse of the light diquark propagator, in momentum 
space the polarization tensor of the scalar diquark reads: 



tr [75 + rrii) 75 (jf^ + m^)] 

4 Jp {{k — pY — mf + ie) {p^ — mj + ie) 

^ {l\m1) + I\m]) + {e - (m, - m,f) dx l\ml{e))^ . (B.6) 



The calculation for the axial vector diquark with different fermion masses m^, mj in the loop can 
be divided into a purely transversal part and an additional contribution proportional {mi—mj)g^^. 
The latter reflects the explicit breaking of 5'f/(3)i? for different quark masses. 



V - « / (27r)4^ 



^ — rrij + ie ^ — ^ — rrii + ie 



{g^uk^-k^k,)7T',+g^,A' . (B.7) 
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Note that in order to respect the Ward identities for the SU{3)f symmetry (as long as = rrij), 
one has to identify the gauge invariant part of vr^^y. The proper-time regularization of the real 
part of the light quark determinant automatically preserves gauge invariance. We obtain: 

= -r^, f dxx{l-x)l\ml{e)) (B.8) 

A'^' = ' ^ ' / dx {{l-x)m,-xmi) l\mUe)) (B.9) 
2 Jo 

with the abbreviation 

m2^.(fc2) ;= xm- + (1 - x)m'^j - x{l - x)k'^ . (B.IO) 

The masses for the light axial vector diquark are presented in Table in the text. 

Note that the polarization tensor for the heavy diquarks is very similar to those for heavy 
mesons 



In momentum space we obtain (see (|T 

(B.ll^ 



The heavy diquark masses following from this can also be found in Table |l| in the text. 
B.2 Baryonic self— energy integrals 

The baryonic self-energy integrals which were used for estimating heavy baryon masses in (^2[) are 
defined as following: 



1 f P P 

Fij{v-p) := / Dij{- - q) g^{- + q) 



H Jq 



Jg {p — qY — {Mf^Y + ie v-q + ie 



^llAv-p) , (B.12) 



Ff^'^ivp) := -- I DUl-q)g,{l + q) 

Z^i v-q + rrij 



H Jq v-p — v-q — AMi + ie q^ — mj + ie 

= hi + _ AM,) Il^{vp - AM,) + ^I' (B.13) 
/i ji 

Ft'^%-p} := -- fD-/--q)g^iP- + q) 

fl J q Z Z 

I f 1 (, 1 q'-{vqr\ 



Jq q^ - {M^jY + ie v-p- v-q + ie \ 3 {M^jY J 
(ll,.{vp) + ll,.{vp)) . (B.14) 



_2l 



18 



References 

[1] M. Voloshin and M. Shifman, Sov. J. Nucl. Phys. 45 (1987) 292, Sov. J. Nucl. Phys. 47 (1988) 
511; 

N. Isgur and M. Wise, Phys. Lett. B232 (1989) 113, Phys. Lett. B237 (1990) 527; 
E. Eichten and B. Hill, Phys. Lett. B234 (1990) 511; 

B. Grinstein, Nucl. Phys. B339 (1990) 253; 
H. Georgi, Phys. Lett. B240 (1990) 447; 

A. Falk, H. Georgi, B. Grinstein and M. Wise, Nucl. Phys. B343 (1990) 1. 

[2] N. Isgur and M. Wise, contribution to Heavy Flavors, A. Buras and M. Lindner (eds.) (World 
Scientific, Singapore, 1992); 
M. Neubert, Phys. Kept. 245 (1994) 259. 

T. Manuel, W. Roberts and Z. Ryzak, Nucl. Phys. B368 (1992) 204. 
J.G. Korner and G. Thompson, Phys. Lett. B264 (1991) 185. 
T. Manuel and G.A. Schuler, Phys. Lett. B349 (1995) 181. 

J.W. Bos, D. Chang, S.C. Lee, Y.C. Lin and H.H. Shih, Phys. Rev. D51 (1995) 6308. 
G. Burdman and J. F. Donoghue, Phys. Lett. B280 (1992) 287. 
M. B. Wise, Phys. Rev. D45 (1992) R2188. 
R. Casalbuoni et al., |hep-ph/960534^ (1996). 

T.-M. Yan, H.-Y. Cheng, C.-Y Cheung, G.-L. Lin, Y. C. Lin and H.-L. Yu, Phys. Rev. D46 
(1992) 1148. 

P. Cho, Nucl. Phys. B396 (1993) 183, ibid. B421 (1994) 683. 

D. Ebert, T. Feldmann, R. Friedrich and H. Reinhardt, Nucl. Phys B434 (1995) 619. 
M. A. Novak, M. Rho and I. Zahed, Phys. Rev. D48 (1993) 4370. 
W. A. Bardeen, and C. T. Hill, Phys. Rev. D49 (1993) 409. 

E. Witten, Nucl. Phys. B145 (1978) 110. 

L. Gamberg, H. Weigel, U. Ziickert and H. Reinhardt, to appear in Phys. Rev. D54 (1996). 

C. G. Callan, I. Klebanov, Nucl.Phys. B262 (1985.) 365. 

J. Schechter, A. Subbaraman, S. Vaidya, H. Weigel, Nucl.Phys. A590 (1995) 655. 

R. Dashen, E. Jenkins and A.V. Manohar, Phys.Rev. D49 (1994) 4713, ibid. D51 (1995) 
2489. 

M.A. Luty and J. March-Russell, Nucl.Phys. B426 (1994) 71. 
R. T. Cahill, C. D. Roberts and J. Praschifka, Aust. J. Phys. 42 (1989) 129. 



19 



[22] H. Reinhardt, Phys. Lett. B244 (1990) 316. 

[23] D. Ebert and L. Kaschluhn Phys. Lett. B297 (1992) 367. 

[24] D. Ebert, H. Reinhardt and M.K. Volkov, in A. Faessler (ed.): Prog. Part. NucL Phys. 33 
(1994) 1; 

R. Alkofer and H. Reinhardt, in Chiral Quark Dynamics, Springer lecture notes in physics, 
1995; 

U. Ziickert, R. Alkofer , H. Weigel and H. Reinhardt, Phys. Lett. B362 (1995) 1. 



[25; 

[26; 
[27; 

[28 

[29 
[30 
[31 
[32; 
[33 
[34 
[35; 
[36 

[37 

[38 
[39 

[4o; 



A. Falk, Nucl. Phys. B378 (1992) 79. 

T. Mannel, W. Roberts and Z. Ryzak , Nucl. Phys. B355 (1991) 38. 

J.G. Korner, M. Kramer and D. Pirjol, Prog. Part. Nucl. Phys. 33 (1994) 787. 

D. Ebert and H. Reinhardt, Nucl. Phys. B271 (1986) 188; 

see also: D. Ebert and M.K. Volkov, Yad. Fiz. 36 (1982) 1265; Z. Phys. C16 (1983) 205. 
N. G. Van Kampen, Phys. Rep. C24 (1976) 171; Physica 74 (1974) 215,239. 
A. Buck, R. Alkofer and H. Reinhardt, Phys. Lett. B286 (1992) 29. 
A. Buck and H. Reinhardt, Phys. Lett. B356 (1995), 168. 
C. Hanhart and S. Krewald, Phys. Lett. B344 (1995) 55. 

N. Ishii, W. Bentz and K. Yazaki, Phys. Lett. B301 (1993) 165; ibid. B318 (1993) 26. 
R.M. Barnett et al. (Particle Data Group), Phys. Rev. D54, 1 (1996) 
CLEO Collab. (P. Avery et al.), Phys. Rev. Lett. 75 (1995) 4364. 

DELPHI Collab., DELPHL95-107, contribution eps0565 to the EPS-HEP95 conference at 
Brussels, 1995; M. Feindt, CERN-PPE-95-139. 



E. Jenkins, |hep-ph/ 9603449 



A.F. Falk, |hep-ph/9603389 



D. Ebert, T. Feldmann, C. Kettner and H. Reinhardt, Z. Phys. C71 (1996) 329. 
H. Reinhardt, Nucl. Phys. A 503 (1980) 825.; 

H. Weigel, R. Alkofer and H. Reinhardt, Nucl. Phys. B 387 (1992) 638. 



20 



